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ABSTRACT 

Multiple  unattended  ground  sensors  are  deployed  for  surveillance,  monitoring  the  movement  of  troops,  military  vehicles,  and 
targeting.  Usually,  the  probability  of  detection  (P D)  of  an  UGS  is  low  and  the  false  alarm  density  (FAD)  is  high.  The 
particle  filter  (PF)  and  range-parametrized  extended  Kalman  filter  (RPEKF)  have  been  used  previously  to  produce  improved 
results  for  the  single  sensor  single  target  bearing-only  tracking  problem  in  the  absence  of  clutter  and  with  unity  probability  of 
detection.  Although  a  great  deal  of  work  has  been  done  for  the  single  target  single  sensor  bearing-only  tracking  problem,  not 
much  is  known  for  the  single  target  multi-sensor  bearing-only  problem  with  low  P D  and  high  FAD.  We  present  algorithms 
and  numerical  results  using  particle  filter  (PF)  and  probabilistic  data  association  (PDA)  for  the  single  target  multi-sensor 
bearing-only  tracking  problem  with  high  false  alarms  and  low  probability  of  detection.  Our  numerical  results  show  that  the 
algorithm  estimates  the  target  state  in  a  robust  manner  in  realistic  harsh  conditions  when  the  probability  of  detection  is  low 
and  the  false  alarm  density  is  high. 

Keywords:  Multi-sensor  Single-target  Tracking,  Unattended  Ground  Sensors  (UGS),  Acoustic  Sensors,  Particle  Filter  (PF), 
Probabilistic  Data  Association  (PDA) 


1.  INTRODUCTION 


Unattended  ground  sensors  (UGSs)  [8], [9], [36]  are  deployed  on  the  ground  and  are  widely  used  in  military  and  industrial 
applications  to  collect  signals  from  remote  sources.  In  the  military  domain,  UGSs  are  used  to  detect,  track,  and  classify 
ground  vehicles  and  aircrafts.  In  the  industrial  domain,  UGSs  are  used  in  mining  operation  that  involves  deep  blasting. 
Common  types  of  UGSs  are  acoustic,  seismic,  seismic/acoustic,  magnetic,  seismic/magnetic,  seismic  string,  passive  infrared 
(PIR),  active  infrared  (AIR),  andbreakwire  [36].  Current  UGSs  are  usually  low-cost,  lightweight,  and  small  size  sensors 
compared  to  the  previous  UGSs. 

Aerial  surveillance  and  human  intelligence  were  primarily  used  to  detect  and  track  the  movement  of  ground  forces  until  the 
mid  1960s.  During  the  Vietnam  War,  it  was  realized  that  these  approaches  are  ineffective  to  monitor  the  force  movement  in 
thick  forests.  Moreover,  long  range  standoff  airborne  radar  sensors  are  not  capable  of  providing  continuous  surveillance  data 
due  to  occlusion  of  the  radar  line-of-sight  by  terrain  in  certain  scenarios.  The  UGSs  can  complement  the  airborne  and 
spaceborne  sensors  to  provide  continuos  surveillance  data  of  the  battlespace. 
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In  this  paper,  we  address  the  problem  of  single  target  tracking  in  clutter  using  multiple  acoustic  UGSs.  An  UGS  determines 
the  bearing  of  a  target  by  processing  the  incoming  acoustic  signal.  In  addition  to  bearing,  other  derived  measurements  such  as 
harmonic  frequencies  are  also  available.  We  restrict  our  analysis  to  single  target  tracking  using  bearing-only  measurements 
from  multiple  acoustic  UGSs.  Usually,  the  probability  of  detection  (  P  D  )  for  an  UGS  is  low  and  the  false  alarm  density 
(FAD)  is  high.  Although  a  great  deal  of  work  has  been  done  for  the  single  sensor  bearing-only  tracking  problem,  not  much  is 
known  for  the  single  target  multi-sensor  bearing-only  problem  with  low  P  D  and  high  FAD.  The  particle  filter  (PF) 
[2],[5],[10],[11],[14]-[17],[30],[32]  and  range -parametrized  extended  Kalman  filter  (RPEKF)  in  [2],  [25]  have  been  shown  to 
produce  improved  results  for  the  single  sensor  single  target  bearing-only  tracking  problem  with  unity  P D  and  no  clutter.  The 
difficulty  in  bearing-only  tracking  is  the  low  information  content  of  the  measurement,  in  addition  to  the  nonlinearity  in  the 
measurement  model.  The  state  estimates  are  also  very  sensitive  to  the  relative  geometry  between  the  sensor  and  target.  Often, 
the  estimates  are  unreliable,  i.e.,  they  have  very  high  variance  [l],[13],[21]-[24],[33]-[35].  We  present  algorithms  and 
numerical  results  using  the  PF  for  the  single  target  multi-sensor  bearing-only  problem  in  clutter  that  represents  realistic  harsh 
conditions. 

In  Sections  2  and  3,  we  present  the  target  kinematic  model  in  two  dimensions  and  acoustic  measurement  model,  respectively. 
We  describe  the  Bayesian  approach  that  provides  a  rigorous  framework  for  state  estimation  with  general  probability  density 
functions  for  the  state  and  measurement  in  Section  4.  Analytic  solutions  using  the  Bayesian  approach  are  not  possible  for 
cases  that  involve  nonlinear  dynamics  for  the  state,  nonlinear  measurement  models,  and  non-Gaussian  distributions.  The 
well-  known  Kalman  filter  estimator  is  a  Bayesian  estimator  for  linear  dynamic  model  with  additive  Gaussian  noise,  linear 
measurement  model  with  additive  Gaussian  noise,  and  prior  Gaussian  distribution  for  the  state.  The  bearing  measurement 
model  is  a  nonlinear  function  of  the  target  state.  The  commonly  used  extended  Kalman  filter  (EKF)  [3], [4], [5], [6], [19]  is  not 
an  optimal  estimator  due  to  the  nonlinear  models  for  the  measurement  and  state  dynamics.  When  the  degree  of  nonlinearity 
is  small  and  the  errors  in  the  measurement  and  state  dynamics  are  small,  the  EKF  is  a  reasonable  approximate  estimator.  The 
EKF  has  been  used  within  other  estimator  configurations  like  the  Interacting  Multiple  Model  (IMM)  estimator  and  the  IMM 
probabilistic  data  association  (PDA)  estimator  to  solve  different  tracking  problems  [3], [4], [6], [7], [27], [28], [31].  When  false 
alarms  are  present  and  the  P D  is  less  than  unity,  we  can  have  zero  or  multiple  measurements  from  a  single  sensor  in  a  scan 
for  a  single  target.  Commonly  used  approaches  for  this  problem  are  the  multiple  hypothesis  tracking  (MHT)  [6], [26]  and 
probabilistic  data  association  (PDA).  The  PF  algorithm  using  the  PDA  [4],[1 1] ,[  1 8]  is  described  in  Section  5.  We  present 
numerical  simulation  and  results  in  Section  6  and  conclusions  in  Section  7. 


2.  TARGET  KINEMATIC  MODEL 


The  nearly  constant  velocity  model  (NCVM)  in  two  dimensions  is  a  widely  used  kinematic  model  [3], [4], [6]  for  ground 
moving  targets.  The  continuous-time  target  state  for  the  NCVM  is  defined  by 


(2.1) 


x(t)  :=  [x(t)  y(t)  x(t) 


where  ( x(t),  y(t))  and  ( x(t),  y(t ) )  are  the  Cartesian  position  and  velocity  coordinates,  of  the  target  at  time  t,  respectively. 
The  continuous-time  NCVM  is  described  by  the  linear  stochastic  differential  equation 

(2.2)  x(f)  =  F  (t)x(t)  +  w(f). 


where 


F(0 


0  0  10 
0  0  0  1 
0  0  0  0 
0  0  0  0 


(2.3) 


(2.4) 


w(f)=[o  0  wx(t)  wy(t)  \ 


F  is  known  as  the  system  matrix  and  wx  (.)  and  wy  (.)  are  mutually  independent,  zero-mean  white  Gaussian  acceleration 


processes: 

(2.5) 

E{wx(t)}  =  0, 

E{  wx  (t)  wx  (x ) }  =  q  x5  (t  -x ), 

(2.6) 

E{Wy(t)=  0, 

E{  wy  (t )  (x ) }  =  q  y  5  (t  - x ), 

(2.7) 

E{wx(t)wy  (r)} 

=  0,  for  all  t,T . 

qx  and  qy  represent  the  power  spectral  density  (PSD)  of  w  Q  and  wy  (•) ,  respectively  [19]. 

We  assume  the  initial  state  of  the  target  to  have  the  prior  distribution 
(2.8)  x(0)  ~  /v(0,P()), 

where  P0  is  the  prior  covariance  of  the  state. 

Sensor  observations  are  available  at  discrete  times  tk,k  =  0,1,... .  In  the  current  work,  we  do  not  address  the  out-of-sequence 
measurement  (OOSM)  that  can  arise  in  a  multi-sensor  tracking  scenario.  Since  the  measurements  are  available  at  discrete 
times,  it  is  necessary  to  discretize  the  continuous-time  model  in  (2.2)  corresponding  to  the  measurement  times.  Discretization 
of  the  continuous-time  dynamics  based  on  this  sequence  of  times  yields  the  following: 

(2.8)  xk+l=<t>kxk+wk. 


where 


(2.9) 


(2.10) 
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yk 

y(fk) 

xk 
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(2.11)  w(  :=  J«I>(tjt+1,?)w(f)£/f,  k  =  0,1,... 

h 

We  can  prove  that  {w  k ,  k  >  0}  is  a  zero-mean  white  Gaussian  process: 


(2.12) 


£{W(t}=0,  E[wkw',]  =  5klQk ,  k,l=  0,1,... 


(2.13) 


wk~N(0,Qk),  k  =  0,1,... 


(2.14) 
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jqx(Atk)3  0  ±qx(Atk)2 
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Qk  is  known  as  the  process  noise  covariance  matrix. 


3.  ACOUSTIC  UGS  MEASUREMENT  AND  CLUTTER  MODEL 


The  number  of  measurements  in  a  scan  can  be  zero  or  greater  than  zero  due  to  PD  <  1  and  the  presence  of  false  measurements 
or  clutter.  Let  mk  denote  the  number  of  measurements  in  the  k'h  scan.  The  set  of  measurements  zk  in  the  kth  scan  is 
defined  by 

(3.1)  **:={*„}£, 

where  z..  is  the  ;th  measurement  in  the  kth  scan.  For  the  acoustic  sensor  z,.  is  a  scalar.  For  the  sake  of  generality,  we 
consider  z ..  as  a  vector.  All  the  measurements  up  to  and  including  the  k th  scan  are  defined  by 

(3.2)  Zk  :={Zj}kj=l. 

We  describe  the  acoustic  sensor  measurement  model,  false  measurement  model,  and  likelihood  functions  for  both  the  target 
originated  measurement  and  false  alarm  in  Sections  3.1,  3.2,  and  3.3,  respectively. 


3.1  Acoustic  Sensor  Measurement  Model 

Acoustic  sensors  are  omni-directional  sensors,  which  determine  bearing  by  processing  the  incoming  acoustic  signal.  The 
bearing  is  a  nonlinear  function  h  of  the  target  state  xk  and  the  sensor  state  sk  at  time  t ,  .  We  assume  that  the  acoustic  sensor 

can  detect  a  target  with  probability  PD  when  the  distance  of  the  target  from  the  sensor  is  less  than  a  maximum  range,  rmax  . 
We  assume  that  the  measurement  noise  at  time  t ,  is  an  additive  independent  Gaussian  noise  process.  Thus,  the  nonlinear 
measurement  model  is  described  by 

(3.3)  zk=h(xk,  sk)+vk,  i(xk-sxk)2 +(yk -syk)2]112  <r 

(3.4)  sk  :=s(4)  =  [sx(4)  sy(tk )  sx(tk)  sy(f*)]  , 

xk  syk  Sxk  S_v*]:=  [sx(li-)  sy(tk)  Lv(ft-)  ^y(^t)]> 


(3.5) 


(3.6) 


E{vk}  =  0,  E{vkvl}=bkIRk,  vk~N(0,Rk), 


(3.7)  Rk=<*L 

where  o  ka  is  the  variance  of  the  bearing  measurement  noise  in  the  k th  scan.  In  practice,  o  ka  can  vary  with  scan  depending 
on  the  distance  of  the  target  from  the  sensor.  For  simplicity,  we  assume  that  the  variance  of  the  bearing  measurement  noise  is 
a  constant  for  all  scans: 

(3.8)  Rk  =a2a. 

We  also  assume  that  the  measurement  noise  V ,  and  the  process  noise  w;  are  uncorrelated  at  all  times: 

(3.9)  E{vk  w;}  =  0,  /or  all  k,l. 


The  bearing  angle  is  measured  from  the  Y-axis  (usually  chosen  along  the  local  North  direction)  in  the  clockwise  direction 
and  the  measurement  function  is  described  by 


(3.10)  h{xk,sk)=  tan  l{xk -sxk,yk -syk),  h(xk, s*  )e  [0,2ji). 

The  sensors  are  stationary  with  respect  to  the  ground  in  our  application.  Therefore,  sxk  =  syk  =  0,  for  all  k. 

3.2  Clutter  Model 


In  addition  to  actual  measurement  originating  from  a  true  target,  there  may  be  false  measurements  observed  by  a  sensor  in  a 
scan.  We  assume  that  the  number  of  false  measurements  mk  in  the  k th  scan  obeys  the  Poisson  distribution.  Thus  the 

probability  of  observing  mk  false  measurements  in  the  k th  scan  is 


-M 


(3.11) 


PfA  (m k)  =~ 


(/V)” 


where  X  is  the  average  number  of  false  measurements  per  unit  area  of  measurement  space  per  scan  and  V  is  the  volume  of 
the  measurement  space.  For  the  omni-directional  acoustic  sensor 

(3.12)  V  =  2 K. 


We  assume  that  the  false  measurements  are  uniformly  distributed  in  the  measurement  space. 


3.3  Likelihood  Function 


The  likelihood  functions  for  target  originated  measurement  and  false  measurement  are 


(3.13) 


P(zk\xk)  = 


j  N[zk',hk(xk,sk),Rk]=  N[zk(h)-hk(xk,sk);0,Rk],  if  zk  is  target  originated, 
[  1  /  V",  if  zk  is  not  target  originated. 


4.  BAYESIAN  ESTIMATION 


A  general  discrete  time  dynamics  for  the  state  x,  e  9?"  of  a  system  is  described  by  [11],[14],[15],[20],[30] 
(4.1)  x*=f*_1(xt_1,wt_1), 


where  f ,x  :  91"  x9i/’  — >  9?"  is  the  system  evolution  function  and  w  k_x  e  9?p  is  a  white  noise  sequence  (known  as  the 
process  noise)  independent  of  the  past  and  current  states.  We  assume  that  the  probability  density  function  (pdf)  of  wA_j  is 
known.  The  function  f;  ,  may  be  a  linear  or  nonlinear  function  of  xk_1  and  wk_x ,  and  the  pdf  of  lwk_x  may  be  arbitrary. 
We  assume  that  measurements  {zk  e  91'" }  are  available  at  discrete  times  and  a  functional  relationship  between  the 
measurement  zk  and  the  state  xk  is  known  [1 1], [14], [15], [20], [30]: 

(4.2)  zk=hk(xk,vk). 


where  hA  :  91"  X  91'  — >  9? m  is  the  measurement  model  function  and  \k  e  9?"  is  a  white  noise  sequence  (known  as  the 
measurement  noise)  with  known  pdf.  We  assume  that  \k  is  independent  of  the  past  and  current  states  and  process  noise. 
The  measurement  function  hA  may  be  a  linear  or  nonlinear  function  of  x/(  and  \k  ,  and  the  pdf  of  \k  may  be  arbitrary.  Let 
Zk  denote  the  set  of  measurements  {zl,z2,...,  zkj  . 


Our  objective  is  to  compute  the  conditional  density  p(xk  I  Zk  )  of  the  state  xk  given  all  the  measurements  Zk  at  time  k. 
Suppose  p{xk_ j  I  Zk~l)  is  known.  Then  p(xk  I  Zk )  can  be  computed  from  p(xk _x  I  Zk~l )  using  the  prediction  and  the 
measurement  update  steps.  Using  the  prediction  step,  the  prior  pdf  of  the  state  at  time  tk  is  given  by  [11], [14], [15], [20], [30] 


(4.3) 


p(xk  I  Zk)  =  | p(xk  I  xk_x,Zk  X)p{xk 4 
=  Jp(xA  I  xk_x)p{xk_x  \Zk~l 


IZ  k~1)dxk_x 
)dxk_x 


where  p(xk  I  x ,  _x )  is  known  as  the  state  transition  density  and  is  determined  by  the  system  dynamics  model  (4.1).  The 
second  step  on  the  RHS  of  (4.3)  is  obtained  from  the  first  using  the  system  dynamics  model  (4.1).  Then  the  prior  pdf  of  the 
state  p(xk  I  Zk  )  can  be  updated  at  time  k  by  Bayes’  rule  using  the  measurement  z ,  [11], [14], [15], [20], [30]: 


(4.4) 


,  p(zk  \xk,Zk~l)p(xk  IZA~1)_  p(zk  \xk)p(xk  IZA-) 
P  "  p(zk  I  Zi_1)  p( zk\ZkA) 


(4.5) 


p(zk  \Zk  *)  =  J  p(zk  I  xk  )p(xk  I  Zk  l)dxk  . 


p(z .  I  x  . )  is  known  as  the  likelihood  function  and  is  determined  by  the  measurement  model  (4.2).  Equations  (4.3)  and  (4.4) 
represent  a  general  recursive  solution  for  the  state  estimation  problem  in  the  Bayesian  framework.  These  equations  are  valid 
for  any  pdf  of  the  state,  process  noise,  and  measurement  noise  with  general  system  dynamics  and  measurement  models. 
However,  closed  form  solutions  are  not  always  possible.  When  and  are  linear,  and  w{  and  \k  are  additive  Gaussian 
noises  with  known  pdf,  (4.3)  and  (4.4)  give  rise  to  the  well-known  Kalman  filter  (KF)  algorithm  [3], [4], [6], [19].  Since 
for  the  acoustic  measurement  model  is  nonlinear,  the  EKF  is  not  an  optimal  estimator  for  the  problem. 

Given  Zk ,  the  conditional  mean  estimator  [11], [14]  ,[15] 

(4.6)  xklk  :=  E{xk  I  Zk }  =  J xkp(xk  I  Zk)dxk 

represents  the  minimum  variance  or  minimum  mean-square  error  (MMSE)  estimator  [3], [4].  The  covariance  of  the 
conditional  mean  estimate  is  given  by  [1 1], [14], [15] 

(4.7)  Pm  :=  E{(xk  -  xm ){xk  - xm )'  I  Zk }  =  J (xk  - xk]k )(xk  - xkik )’ p(xk  I  Zk  )dxk . 

We  emphasize  that  if  the  posterior  distribution  is  multi-modal,  then  the  conditional  mean  and  associated  covariance  are  not 
sufficient  statistics. 


5.  SINGLE  TARGET  TRACKING  IN  CLUTTER  USING  PARTICLE  FILTER 


The  computation  of  the  likelihood  function  in  the  absence  of  clutter  is  straightforward.  We  use  the  PDA  approach 

[4],[  1 1] ,[1 8]  to  compute  the  likelihood  function  required  by  the  PF.  First  we  present  the  algorithm  to  compute  the  likelihood 

function  and  then  we  present  the  steps  of  the  PF  algorithm. 

5.1  Likelihood  Function 

When  there  are  m k  measurements  in  the  k th  scan,  we  have  the  following  mutually  exclusive  and  exhaustive  hypotheses  [4]: 

f[z  kj  is  the  target  originated  measurement],  j  =  1,2  k , 

(5. 1)  0  kj  =  ] 

[  [none of  the  measurements  is  target  originated],  j  =  0. 

Using  the  total  probability  theorem  [3] 

mk 

(5.2)  p( zk  \xk)  =  ^p{zk  I Qkj,xk)P(Qkj  \  xk), 

i= o 

where  P(Qkj  lxir)  and  P(zk  I  ®kj  >xk)  are  the  probability  and  the  likelihood,  respectively,  of  the  hypothesis  Qkj  .  Since  the 
measurements  are  conditionally  independent, 

mk 

(5.3)  p(zk  \$kj,xk)=Y\p(zkr  ^kj^k)- 

r= 0 


Using  the  measurement  model  (3.3)  and  likelihood  in  (3.13) 


(5.4) 


P(Zkj  I  ®kj’Xk)  = 

j  N[Zy-,hk(xk),Rk]  =  N[zkj(k)-hk(xk);0,Rk],  if  zkj  is  target  originated. 


[l/U, 


if  z  kj  is  not  target  originated. 


where  V  is  the  volume  of  the  measurement  space.  For  generality,  we  treat  zkj  as  a  vector  measurement.  Using  (5.4) 
(5.3),  we  get 


(5.5) 


p{zk  1 9  j(k),xk  )  = 


\v~mt+1N(?kj-,0,Rk), 

[  V~n'k, 


j  =1,2,...,  mk, 

7=0, 


(5.6)  ?kj  :=zkj-hk(xk), 

(5.7)  #(?#;  0,Rk)  =  ^L, 

Vk 

(5.8)  e.W-expf-I?;.^)^-1?,.], 

(5.9)  n,:=[(27t)'"|R,|]l/2. 

For  simplicity  in  notation,  define 

(5.10)  jj:=P(Qj(k)\xk). 


Then  we  get  [4] 
(5.11)  Y  j  = 


pD+a-pD) 


PfA  (mk  ) 


j  =  1,2,..,  mk 


(1  ~Pn) 


P FA  (,nk) 


PFA  (mk  -1) 


pD  +  a-pn) 


PfaW 
PFA(mk  ~  1) 


7  =  0. 


Using  the  Poisson  model  (3.11)  for  clutter,  we  get 

(5.12)  W  _(»>,-l)!aU)m‘ 

pfa  (mk  - 1)  mk !  (A,  V)'"1  _1  "U 

Substituting  (5.12)  in  (5.11),  we  get 


(5.13) 


\pD[mkPD  +(1  -PD)ly]~\  j  =1,2,...,  mk 
1  \(l-PD)Xv[mkPD+(l-PD)Xv]-\  7=0. 


Define 


(5 . 14)  Pit,-  :=  p( z k  10*,-, x,)P(0^  I  x*),  j  =  0,1,...,  mk . 
Then 

(5.15)  p  >-^v(V  o,R*)y,, 

J  l  Yj»  7=0. 


Using  (5.7)  and  (5.13)  in  (5.15),  we  get 


(5.16) 


e/g,  j=l,...,mk, 
»,  1-0. 


(5.17)  bk  :=\(\  -  PD  )/  PD  \  (Kr\k ), 

where  ^  is  a  constant.  Substitution  of  (5.16)  in  (5.2)  gives 


mk  mk 

(5.18)  p(  zk  I  x4)  =  £„(«*  I0a7,xa.  ^(O^  lx*,)  =  £[&* 

7=o  7=i 

5.2  Steps  of  PF  Algorithm 

1.  Initialization 

Set  scan  index  k  =  0. 

Sample  {x^  ~  p0(x)}^i  • 

2.  Increment  scan  index  k:  k=k+ 1. 

3.  Prediction 

Generate  N  samples  of  the  process  noise: 

(5.19)  { w‘(k,k-l)~  N(0,Q(k,k-l)}fLv 
Compute  N  predicted  state  vectors: 

(5.20)  {x-  =<D(fe,fc-l)x^1  +vr‘(k,k-l)f=l. 

4.  Compute  likelihoods  [p( zk  I  x'k )  . 

(5.21)  v kj  :=  z kj  —  h*  (x* ). 

We  note  that  the  predicted  measurement  is  the  same  for  all  measurements  { . 


(5.22) 


4  =exp[-^(?ip/R^14]=exp[-|[zy-h*(x4],Ri1[z^-h-t(xi)]]- 

m. 

(5.23)  /£=[&* +£4-],  i  =  1,2,..., 


,  AT. 


7=i 


5.  Update  weights  and  normalize 

i  J  i 

(5.24)  w',  -  4^'  k  .  i  =1,2,...,  A/. 

Y_,w'k-\lk 

1=1 

6.  Compute  the  measurement  updated  state  estimate: 

N 

(5.25)  =Y^wkxk 

1=1 

7.  Compute  effective  sample  size  (Neff  ) 

N 

(5.26)  AT*  =1  /£(wi)2. 

1=1 

^eff  >  ^thres  -  then 

Go  to  step  2. 

Else 

8.  Resample  a  new  set  { x)  [T,  by  sampling  with  replacement  Atirnes  from  the  discrete  set 

U*}£i  where  Pr(xf=x/)  =  w/. 

Set  {<  =  1/A0m- 
Go  to  step  3. 

End 


6.  NUMERICAL  SIMULATIONS  AND  RESULTS 


The  truth  trajectory  of  the  target  and  the  positions  of  four  acoustic  UGSs  are  shown  in  Figure  6.1.  We  use  the  nearly  constant 
velocity  model  to  generate  the  truth  trajectory.  The  initial  position  and  speed  of  the  target  are  (-220  m,  300m)  and  40  km/hr, 
respectively.  The  initial  azimuth  angle  of  the  target  is  45  degrees.  The  PSD  of  each  component  the  acceleration  process  noise 
( q )  for  the  truth  trajectory  is  0.01  m2s'3.  We  use  5000  particles  throughout  our  simulation.  The  detections  and  false 
measurements  at  various  scans  for  the  case  when  P D  =  0.9  and  the  average  number  of  false  measurements  per  scan  is  one, 
are  shown  in  Figures  6.2.  The  truth  and  PF  estimated  position  and  velocity  for  a  benign  scenario  where  the  P D  =  0.9  and 

average  number  of  false  measurements  per  scan  is  one,  are  shown  in  Figures  6. 3-6. 5.  We  observe  from  Figures  6. 3-6. 5.  that 
the  PF  algorithm  estimates  the  position  and  velocity  of  the  target  accurately  for  this  benign  scenario.  The  root  mean  square 
(RMS)  position  and  velocity  errors  are  presented  in  Table  6.1.  We  present  results  for  other  cases  in  Table  6.1  by  varying  the 
P D  and  average  number  of  false  measurements  per  scan.  We  observe  from  Figures  6. 3-6. 5  and  6. 7-6. 9  that  the  large  RMS 

errors  are  mainly  due  to  transients  in  the  early  part  of  estimation  process.  After  the  transient  phase,  the  PF  estimated  state  is 
close  to  the  true  state. 
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Figure  6. 1.  Sensor  Locations  and  Truth  Trajectory  of  the  Target  using  the  Nearly  Constant  Velocity  Motion 


Figure  6.2.  Detection  and  false  alarms  at  various  scans  for  high  probability  of  detection  and  low  clutter. 


Table  6.1.  PF-estimated  position  and  velocity  RMS  errors  using  120  scans. 


Bearing 

Sigma 

(deg) 

PD 

Average 
Number 
of  False 

Alarms 
per  Scan 

RMS 

Position 

Error  (m) 

RMS 

Velocity 

Error  (m/s) 

3 

0.9 

1 

38.824 

3.279 

3 

0.7 

1 

67.682 

3.315 

3 

0.7 

2 

94.076 

3.654 

3 

0.7 

3 

132.235 

4.118 

Figure  6.3.  Truth  and  PF-estimated  trajectories  for  high  probability  of  detection  and  low  clutter. 


Figure  6.4.  Truth  and  PF-estimated  X-component  of  velocity  for  high  probability  of  detection  and  low  clutter. 


Figure  6.5.  Truth  and  PF-estimated  Y-component  of  velocity  for  high  probability  of  detection  and  low  clutter. 


Figure  6.6.  Detection  and  false  alarms  at  various  scans  for  low  probability  of  detection  and  high  clutter. 


Figure  6.7.  Truth  and  PF-estimated  trajectories  for  low  probability  of  detection  and  high  clutter. 


Figure  6.8.  Truth  and  PF-estimated  X-component  of  velocity  for  low  probability  of  detection  and  high  clutter. 


Figure  6.9.  Truth  and  PF-estimated  Y-component  of  velocity  for  low  probability  of  detection  and  high  clutter. 


7.  SUMMARY  AND  CONCLUSIONS 


In  this  paper,  we  have  addressed  the  single  target  multiple  acoustic  UGS  tracking  in  clutter  using  the  particle  filter  (PF) 
algorithm.  We  have  used  realistic  values  for  the  probability  of  detection  and  false  alarm.  We  have  demonstrated  that  the  PF 
algorithm  works  in  a  robust  manner  when  the  probability  of  detection  is  low  and  the  false  alarm  is  high  as  is  the  case  in 
realistic  harsh  scenarios.  In  our  future  work,  we  plan  to  compare  the  performance  of  the  PF  with  the  EKF  using  the  PDA 
approach  and  analyze  the  estimation  accuracy  by  varying  the  accuracy  of  the  acoustic  sensor  measurement. 
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